Gravitational collapse of a spherical star with heat flow 
as a possible energy mechanism of gamma-ray bursts* 
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We investigate the gravitational collapse of a spherically symmetric, inhomogeneous star, which 
is described by a perfect fluid with heat flow and satisfies the equation of state p = p/3 at its center. 
In the process of the gravitational collapsing, the energy of the whole star is emitted into space. 
And the remaining spacetime is a Minkowski one without a remnant at the end of the process. For 
a star with a solar mass and solar radius, the total energy emitted is at the order of 10 54 erg, and 
the time-scale of the process is about 8s. These are in the typical values for a gamma-ray burst. 
Thus, we suggest the gravitational collapse of a spherical star with heat flow as a possible energy 
mechanism of gamma-ray bursts. 
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More than 30 years ago, Klebesadel, Strong and 
Olson announced the discovery of the gamma-ray 
bursts [1] . Since then, thousands of gamma-ray bursts 
have been observed [2]. fn particular, since the ad- 
vent of the Compton Gamma-Ray observation in 1991, 
the gamma-ray bursts have been detected daily. The 
durations of the gamma-ray bursts range from about 
30ms to over 500s [3] . The inferred isotropic luminosi- 
ties are typically on the order of lfj 51 ~ 52 erg/s [4] . The 
abundance of observations has led to a well described 
gamma-ray burst phenomenology. 

In this Letter, we study the gravitational collapse 
of a spherical star with heat flow in a generalized 
Friedmann coordinate system. It is supposed reason- 
ably that a fluid star is spherically symmetric but in- 
homogeneous and that the equation of state at its cen- 
ter is p = p/3. We find that all energy of the star will 
be emitted in the process of collapse, and an empty 
flat spacetime will be left behind. For a star with 
about a solar mass and a solar radius, the energy at 
the order of 10 54 erg will be emitted into space within 
about 8s. It is the typical value for a gamma-ray burst 
in an isotropic model. Thus, we suggest the spherical 
gravitational collapse with heat flow as a possible en- 
ergy mechanism of gamma-ray bursts, which is a dif- 
ferent mechanism from the ones discussed extensively 
in literature [5]. 

In the pioneer paper to understand the late stages 
of stellar evolutions [6], the collapsing star is sup- 
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posed to consist of homogeneous, spherically symmet- 
ric, prcssurclcss, perfect fluid and to be surrounded 
by an empty space. The interior of the star may be 
described by the Friedmann- Robertson- Walker metric 
[7] 

ds 2 = dt 2 - a 2 (i) ( + r 2 dn 2 ) , (1) 
\ 1 — kr 2 J 

where d£l 2 = d9 2 + sin 2 9dip 2 is the metric on a unit 
2-sphere, k is times of the initial energy density 
of dust in the unit of c = 1, and a(t) is the solution of 
a 2 {t) = kia-^t) - 1]. 

In an astrophysical environment, a star usually 
emits radiation and throws out particles in the pro- 
cess of gravitational collapse. In this situation, the 
heat flow in the interior of a star should not be ignored 
and the exterior spacetime is no longer described by a 
Schwarzschild metric. To take the radiation of a star 
into account, the interior solution of the gravitational 
collapse of radiating stars should match to the exterior 
spacetime described by the Vaidya solution [8] 

dg2 = (1 _ 2MW )(i(;2 + 2dydR _ R 2 dQ2> (2) 
it 

which has been studied extensively [9-11]. In particu- 
lar, the gravitational collapse of a radiating spherical 
star with heat flow has been studied in an isotropic 
coordinate system [11] 

ds 2 =dt 2 - B 2 (t,r)(dr 2 +r 2 dQ 2 ), (3) 

where B{t,r) = b 2 (t)[l - A(i)r 2 ] -1 with suitable b(t) 
and \(t) functions. They called the solution the 
Fricdmann-likc solution. 

We study the gravitational collapse of a spheri- 
cally symmetric, inhomogeneous star in a generalized 
Friedmann coordinate system 

ds2 = df2 - a2 Kr|p +rW ) (4) 
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with k being a function of t. 

The stress-energy tensor of the fluid with heat flow 
but shear-free and without viscosity is given by 

= (p + p)u»u v - pg^ + q^u" + q^u^ , (5) 

where p and p are the proper energy density and pres- 
sure measured by the comoving observers respectively, 
u M is the 4- velocity of the fluid, and q^ is the heat flow, 
and q 11 satisfy 

u^u 1 * = f and <Z^u M = 0. (6) 

For the spherically symmetric collapse, one has 

u" = («°,« 1 ,0,0) and q" = (q°,q 1 , 0,0). (7) 



The Einstein's held equations 



(8) 



and the covariant conservation of stress-energy tensor 
give rise to 
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with the identity 
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where the dots denote the derivatives with respect to 
time t, and 

A = - 1 (16) 
y/l - k(t)r 2 

X = 3aaAA + a 2 AA, (17) 
Y = [2{k + d 2 - ad)A - aaA - a 2 A]A, (18) 

Z = —aA. (19) 
r 

X, Y, Z should satisfy the following solvable condition, 
Z 2 > AXY. (20) 



For spherical collapse, the dominant energy conditions 
[12] require 



which lead to 



P -p > 0, 

P + p>2^/-q^, 



X + Y > \Z\ > 0. 



(21) 
(22) 



(23) 



At the center of the star r = 0, A, A = and thus 
X, Z vanish. Hence, 

Mo|r=0 = 1) «l|r=0 = 0, 9 |r=0, 9i|r=0 = 0, (24) 
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(25) 
(26) 



The exterior solution outside the collapsing star 
is described by Vaidya metric (2). The interior solu- 
tion and exterior solution should satisfy the junction 
conditions on the surface of the star [13], 



R s 
dr 
di 

pl = 



= -a- 2 A- 2 ^ 
u 



(27) 
(28) 
(29) 



and 



[dt 2 -a 2 A 2 dr 2 ] s 

- I dv 2 + 2dvdR 
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The Misner-Sharp mass [14] for the star is then 



M(v) 



2G 



(30) 



(31) 



To make further calculations, the equation of state 
for the matter of a star at r = needs to be assigned. 
Without loss of generality, we choose <z(0) = 1. Af- 
ter the equation of state is assigned, one may obtain 
the expression for a(t), k(t), p(r,t), p(r,t) and r s (t) 
by solving Eqs.(25), (26), (9), (10), (15), (28), (29). 
Then, one may determine the Misner-Sharp mass 
M(v) and the surface radius R s from Eqs.(31) and 
(27). Fig. 1 presents the numerical solution in the case 
p\ r= o = p/3| r= o and the Uq solution in Eq.(ll) with 
the initial conditions r s (t = 0) = R s (t = 0) = Rq, 
d(t = 0) = 0, R s (t = 0) w 0, which means the sur- 
face of star is almost stationary in the initial state. 
Since p and p are determined by Eqs.(9) and (10), 
the initial state of the star is not homogeneous and 
the equation of state in the star is generally deviated 
from the radiation. In the figure, the time and radius 
r s are in the units of Rq/c and Rq, respectively. In 
the process of evolution, r s (t) (dash curve) decreases 
monotonically while a(t) (dotted curve) almost keeps 
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FIG. 1: Numerical solution of the gravitational collapse of 
spherical star with heat flow in the case p| r =o = p/3| r =o 
and the solution in Eq.(ll). The horizontal axis is the 
time t in the unit of Rq / c. In the process of evolution, r s (t) 
(dash curve, in the unit of Ro) decreases monotonically 
while a(t) (dotted curve) keeps almost a constant. 2M/R 3 
(real curve) in the unit of c 2 /G increases first and then goes 
to as R 3 — > 0, which implies that the star disappears at 
the end of the collapse. 

a constant. 2M/R S (real curve) in the unit of c 2 /G 
increases first and then goes to more quickly than 
R s itself as R s — > 0. It implies that the star will 
disappear without the appearance of a horizon in the 
process, leaving Minkowski spacetime behind. 

Fig. 2 shows the evolution of p and p at the bound- 
ary, in unit of c 4 /(8itGRq). The emission of the star 
arrives its maximum value at about 0.9Rq/c after the 
beginning of the collapse. In the late stage of the pro- 
cess, the equation of state at the boundary as well as 
at each point in the star tends to p — At the end 
of the process both the energy density and the pres- 
sure become 0, which confirms that the whole star is 
radiated out into space in the process. 

Furthermore, our numerical analysis shows that 
the star with p| r= o = -|p|r=o will radiate its whole 
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FIG. 2: Evolution of p and p at the boundary. The hori- 
zontal axis is the time in the unit of Ro/c. p and p axe in 
the unit of c 4 /(8tyGRo). 



mass in the process of the collapse, without the ap- 
pearance of a horizon, for different initial values. 

In conclusion, we have presented a solution of 
the Einstein equations, which describes the gravita- 
tional collapse of an inhomogeneously-distributed but 
spherically-symmetric star consisting of perfect fluid 
with heat flow but without viscosity. Wc have shown 
that if the star satisfies Eqs.(9) and (10) and p = p/3 
at the center, the whole star will be emitted in the 
process of the gravitational collapse in about 3.3i?o/c 
interval, without the formation of a horizon, and a 
Minkowski spacetime is left at the end of the process. 
For a star with about a solar mass Mq and a solar ra- 
dius Rq, a huge mount of energy (about 1.8 x 10 54 erg) 
will be emitted into space within 7.66 seconds, which 
are typical values for a gamma-ray burst. Thus, the 
gravitational collapse of such a star might provide an 
energy mechanism for gamma-ray bursts. 

Since heat is flowing across the boundary of the 
collapsing sphere and the pressure of the fluid at the 
boundary is not zero, the exterior space cannot be con- 
sidered as vacuum. We have adopted a Vaidya space 
as the exterior of the spherical collapsing stars. In the 
Vaidya space, massless radiation is dominated. The 
spectra of the radiation from collapsing stars depend 
on the equation of state of interior matter. We wish 
the combination of the model with phenomcnological 
studies on possible evolution process of the equation 
of matter state can give a good fit to the spectra of 
the gamma-ray bursts. 

As a solution of Einstein equations, the advantage 
of the model is that it transmits all initial energy of 
a star into radiation and then supplies up to 10 54 erg 
radiated energy by a popular star with solar mass. 
In this collapsing scenario, the gamma-ray bursts are 
supposed to be spherical. Thus, one need not try to 
form beaming of energy, whose dynamical mechanism 
is not well understood up to now. In fact, the most 
important motivation of the study is try to find a pos- 
sible candidate for the replacement of the models for 
the gamma-ray bursts including beaming. Of course, 
other ways should not be ignored. 

In the model, the initial energy, which is trans- 
formed into radiation and emitted completely in the 
process, and initial radius of a star are only param- 
eters of the solutions. Different initial energy and 
radius will result in different energy release and dif- 
ferent collapsing timescale. Therefore, the model for 
the energy mechanism of gamma-ray burtsts will sup- 
ply a wide spread in total energy release and various 
timcscalcs. 

Finally, it should be stressed that if the solution 
really serves as a model for gamma-ray bursts, fur- 
ther studies on many problems are needed, which have 
been out of the coverage of the letter, such as, the 
trigger of collapse, the mechanism of energy transfor- 
mation, the constraints on the spread in total energy 
release and timescale in the observation, light curve, 
spectrum of emmission, and so on. 
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